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Q-i! Abstract 



We present the exact 0(a) correction to the process e + e~ — > 
e + e~ + 7 in the low angle luminosity regime at SLC/LEP energies. 
We give explicit formulas for the completely differential cross section. 
As an important application, we illustrate the size of the respective 
corrections of 0(a 2 ) to the SLC/LEP luminosity cross section. We 
show explicitly that our results have the correct infrared limit, as a 
cross-check. Some comments are made about the implementation of 
our results in the framework of a Monte Carlo event generator. This 
latter implementation will appear elsewhere. 
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1 Introduction 



Recently, new luminometers at LEPjIj] have made measurements of the lu- 
minosity process e + e~ — > e + e~ + 71(7) at the experimental precision tags 
below .1%. This should be compared with the prediction by two of us (S.J. 
and B.F.L.W.), with E. Richter-Was and Z. Was, of these processes at the 
.25% precision tag in Ref. using the YFS Monte Carlo event genera- 
tor BHLUMI2.00. Recently, using version 4.0 of BHLUMI, the authors in 
Ref. J2| have reported precision tags ~ .1% for the ALEPH SICAL detector's 
asymmetric acceptance theoretical cross section prediction || |J and .16% 
precision for the general ALEPH SICAL acceptance. The four of us, with 
W. Placzek, E. Richter-Was, M. Skrzypek and Z. Was 0, have reported 
the entirely equivalent conservative theoretical precision tag ~ .15% for the 
bremsstrahlung correction for the general acceptance of the ALEPH SICAL 
detector. These results are currently being extended to the other new LEP 
luminometers 0. It is clear that the precision tag on the theoretical predic- 
tion of <T£, the SLC/LEP luminosity cross section, urgently needs further 
improvement over the ~ .15% level if the theoretical uncertainties are going 
to be reduced to the required level of one half of A<j£ Xp , the respective ex- 
perimental uncertainty, so that they do not obscure the comparison between 
theory and experiment in the high precision Z° physics tests of the Standard 
Model of the electroweak interaction. 

From Table 2 in Ref. and from Table 1 in Ref. ||, we see that the 
subleading part of the 0(a 2 ) pure bremsstrahlung correction to remains, 
at this writing, a dominant part of the outstanding theoretical uncertainty, 
where it contributes directly to the physical precision and indirectly to the 
technical precision both at a level ~ .1% itself. Exact results on the double 
bremsstrahlung process itself have been derived in Refs. |7j and || and, in 
fact, the entire leading log part of the 0(a 2 ) correction to o~£ has recently 
|9| been incorporated into BHLUMI4.0. What remains to be done rigorously 
then is to compute the remaining exact subleading part of the bremsstrahlung 
correction to ac and incorporate it into our YFS Monte Carlo event generator 
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BHLUMI4.xx to check precisely the size of this effect in &£. 

In this paper, we present the exact results which are necessary to evaluate 
this remaining unquantified sub-leading part of the bremsstrahlung correc- 
tion to o~£ with particular emphasis on the acceptances of the new LEP 
high precision luminometers. While partial results and estimates on this 



sub-leading correction to ac have appeared elsewhere fTO], [11]] , our results are 
the first ever, fully differential exact results of their type. Their detailed 
implementation into BHLUMI4.xx will appear elsewhere [[]]. 

More precisely, the lone outstanding contribution from bremsstrahlung 
processes to which dominates the physical precision part of the error in 
BHLUMI4.0 from 0(a 2 ) bremsstrahlung is the sub-leading part of the virtual 
correction to the single bremsstrahlung process. Thus, it is this latter process 
which we shall compute exactly in what follows in a fully differential manner, 
as is needed for Monte Carlo event generator applications. We repeat — such 
a completely differential, exact 0(a 2 ) single bremsstrahlung calculation has 



not appeared elsewhere. (See, however, Refs. |nj [TT| for various levels of 
partial results.) 

Our work is organized as follows. In Sect. 0, we set our kinematic and 
notational conventions. In Sect. ^, we analyze the processes of interest to 
us using the algebraic program FORM JL^]. In Sect. ||, we present numeri- 
cal results which illustrate checks on our work in the SLC/LEP luminosity 
regime. Sect. |5] contains our summary remarks. The Appendices contain 
some technical details. 



2 Preliminaries 

In this section we set our kinematical notation and calculational conventions. 
We begin with the kinematics. 

The process under discussion is illustrated in Fig. [I] We consider the 
one-loop corrections to the process e + (pi) + e~(qi) — > e + (p 2 ) +e~(g 2 ) +j{k) 
in the low angle regime of the SLC/LEP luminometers, where 8 e +,8 e - G 
[25 mrad, 70 mrad] if 9 e +,9 e - are the CMS scattering angles of e + , e~ in 
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the Z° resonance energy regime. The kinematics is illustrated in Fig. |T|. It 



can be shown [|13| , |T0|| that, in this low angle regime, graphs involving the 
exchange of more than one virtual photon line between different fermion 
lines are suppressed in the 0(a 2 ) correction to the cross section. Thus, we 
do not need to calculate these graphs for the 0(a 2 ) corrections of interest 
to us here. Further, it can be shown that those terms in the cross section 
involving interference of the radiation from the e + line with that from the 
e~ line, in the low-angle Bhabha scattering (LABH) regime, are severely 
suppressed as well — this is the so-called up-down interference suppression 
fl"3| . Hence, we will only need to calculate the graphs, as those shown in Fig. 
[I], where only one photon is exchanged between e + and e~ lines and in which 
a virtual correction exists on one of these lines. There are a total of 36 such 
graphs, excluding vacuum polarization, and the ten electron-line emission 
graphs giving non-trivial results in our on-shell renormalization scheme are 
shown. The associated positron-line emission graphs must be calculated as 
well. The s-channel exchange graphs will also be calculated for completeness 
and we will see that they are indeed negligible at the level of accuracy of 
interest to us here. 

For the actual calculation, we rely on two technical tools. We evaluate 
the general 7 emission amplitudes using the formulation of the CALKUL 
Hp , |15fl methods given by Xu et al. fllBf . Thus, our massless fermion spinors 
\p, ±), of helicity ±, are such that their spinor product is 



(p\p'h = &-\p\+) = (p x +ip y )^-(p' x + ip' y )^ (2.1) 

VP+ yjp' + 

W)- = (p'\pT+ 

for p + = p° + p z in an obvious cartesian coordinate notation for the 3- 
momentum p (i.e. p^ = (p° , p 1 , p 2 , p 3 ) = (p°,p) = (p°,Px,P y ,Pz)), and our 
corresponding photon polarization vectors are, for helicity p, 
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Figure 1: 0(a 2 ) single bremsstrahlung correction in e + e — > e + e at low 
angles. Only electron line emission graphs are shown. 
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where k is the photon momentum and h is an auxiliary massless 4- vector. As 
usual, h may be chosen to simplify a given gauge-invariant set of graphs flT6|. 



For a given helicity, three of the ten graphs in Fig. [1] can be eliminated in 
this manner, leaving seven to be evaluated. The methods of Xu et al. are 
augmented for the evaluation of virtual corrections by the algebraic manip- 
ulation program FORM [|i~2| , which we use to evaluate our typical one-loop 
integrals via reduction to the now-standard scalar integrals; here we use the 
realization of the standard scalar integrals of Refs. |T7|, |I8||, since they were 



numerically stable enough for our applications. It was with the aid of these 
two computational techniques that we have calculated our results for the 
0(a) correction to the single bremsstrahlung process e + e~ — > e + e~ + 7 at 
low SLC/LEP luminosity regime angles. Our results are presented in the 
next section. 

3 Exact 0(a 2 ) Results for e + e" — > e + e~ + 7 

In this section we present our results for the 0(a) correction to the process 
e + e~ — > e + e~ + 7 at low angles. We discuss first the 27 bremsstrahlung 
effect. 

For our 27 bremsstrahlung correction, we use the results of three of us 
(S.J., B.F.L.W, and S.A.Y.) in Ref. ||, where the exact result for e + e~ — > 
e + e" + 27 was computed with the methods of Xu et al. JL6]. These results 
have been checked in Ref. [|TU . 



For the virtual correction to e + e~ — > e + e~ + 7, we organize our results 
in terms of the amplitude for 7 emission from the electron line and for emis- 
sion from the positron line, neglecting the so-called up-down interference 



terms [jT3| , since these are known to contribute negligibly at the level of our 
current precision of interest. 

The electron-line emission amplitude with one virtual photon is 

A e w = (Wo + ^1X1 + F2I2) (3.3) 
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where 



2V2p 



ipi,P2)-p({qj,Pi) P ) 
(qi,k)p(q2,k) p 



(3.4) 



X 



2y/2fi 



2^2p 



(gj, fc)-p(P2,gj)-A(gj,Pl)A 
(gj, fc)p(gl,g2>-p 

(q h k)- p {p 2 ,qj)-\{qj,pi)x 
(gj, k) p (qi,q 2 )-p 



(3.5) 



(3.6) 



where the helicity-dependent indices i, j, J are given by 



^ = {2 ifp = ±A, (j,J) = { 



(1,2) 
(2,1) 



if p = ±p, 



(3.7) 



where A, p, p are the helicities of the positron, electron, and photon. The 
function X is proportional to the electron line Born amplitude: 



Fo(p = p) = -8 - 8m 2 e C^ 3 + ri (t q - n)- 1 - 2{t q C 12A + t p Cfl) 

- {t q C 124 -nCti -{t q + r 2 )C\il +( ri -r 2 )C 2 34 + ^[234} 

x tqr^in -r 2 )(t q -n)- 1 

+ {t q C l2i + r 2 C^ r 3 2] -(t q -n)C[Z 2] +(ri-r 2 )C 2U - r 2 t q D\^ } 
+ QB 12 + (B [ £ ] -B 34 )r 1 (t q -r 1 )- 1 {l-3t q (t q + r 2 )- 1 } 

- 6^34 + (B 24 -B u ){2t q r 1 (r 1 -r 2 )' 1 (t q -r 1 )- 1 }, 



Born 




(3.8) 



The form factors are 



(3.9) 
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Tx{p = n) = 2t g (n - r 2 )- x - t q (t q - n)- 1 
+ {t q C 1M -nCfe 1 -{t q + r 2 )Cti\ + (n-r 2 )C 2U + r^D&L } 

x{t p t q r 2 2 (t q - ri) _1 (t g - r 2 ) + §<J P)1 } 
- {^C 124 + r 2 ci 23 r21 - (*, - ri)Ct? ] + (n - r 2 )C 234 - ^[234 } 

+ (B [ ; 3 l] - B 3i ) t p t q {t q - n)" 1 {2r 2 - x - 3(t, + r,)- 1 } 

+ 2(E 24 - B u ) t p t q {n - r 2 )- 1 {(r 1 - r^ 1 - t^J - n)" 1 }, 

(3.10) 



T 2 {p = fi) = -2t q { ri - r-2)- 1 + + r,)- 1 
- {t q C l2A -nCti -(t q + r 2 )Ct£ + (r 1 -r 2 )C 234 +nt q D^i} 

x{t p t q r2 2 + -nr^dp^} 
+ {t q C 1M + r 2 C [ ^ ] - (t q - ri )C l 13 r A 2] + (r 1 -r 2 )C 234 - r^D^ 1 } 

x {R»,-i} 

+ 2(5 24 - 5 34 ) t p t q (n - r,)- 1 {r 2 - x - (n - r,)- 1 }. 

(3.11) 

The opposite helicity cases may be obtained from the above results using 
the substitutions (for i — 0, 1, 2) 

J"i(p = -n, n, r 2 ) = Ti{p = /J,, -r 2 , -n) . (3.12) 

The scalar integrals -B, C, -D are defined in Refs. ]T7| (we evaluate them 
using algorithms from Refs. [|l7j]) and in Appendix [A]. The kinematic vari- 
ables s, s , tp, t q , Ti are defined as 

s = (pi + qi) 2 , s' = (p 2 + q 2 ) 2 , 
t p = (P1-P2) 2 , t q = {qi~q 2 ) 2 } (3.13) 
Ti = 2qi ■ k. 
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We have found that at low angles in the SLC/LEP energy regime, the 
JF X terms in ( |3.3| ) are often a good approximation to the entire result, and 
that these terms are in turn well approximated by the simple expression 

i.p 

A e~ _ ^ T approx T ("JItO 

A (l) approx - 167T 2 tp ^0 M , l^ i4 J 

^approx = _ 8 _ 8m 2 C m 3 _ 2( j. qCi ^ +tpC m 4)+6Bi2 _ 6BM 

3 mi mi mi 

+ 41n^4ln^ + 41n^4ln^-81n^. 

■mi m e mi m e m e 

(3.15) 

This approximation has been compared to the complete result in detail and 
we have found it to be within 10% of the result fl3~3] ) throughout most of the 
final particle phase space; such comparisons will be presented in more detail 
elsewhere H. 



The analog of ( |3.3j ) for positron line emission is obtained by crossing: 

A (i) = A (i) (Pi ^ -?2>9i -P2, A -fi). (3.16) 
The differential cross section associated with (13.31) is the usual 



dnldkdn k ~ Wyb J p Re + A ^ ^ l+orn + ALrn ) ' (3 ' 17) 

where Q is the outgoing positron solid angle in the lab frame, flk is the photon 
solid angle, and the up-down interference terms can be neglected here. 

We have used massless fermion spinors to calculate Q3.31 ) - ( |3.17|) . The rel- 



evant mass effects can be restored by the standard methods already published 
in Ref. [ I5| . We have done that, and it amounts to adding an m e 2 -correction 



term for each external emission line in Fig. [l|. The corresponding corrections 
are given in the Appendix || for completeness. In this way, we have arrived 
at an exact 0(a) correction to the process e + e~ — > e + e~ + 7 with all relevant 
mass effects taken explicitly into account. 
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4 Results and Checks 



Our results ( ^.3| ) - ( |3.17| ) are readily introduced into the Monte Carlo pro- 



gram BHLUMI 0, [I of two of us. This will be presented elsewhere Here, 
we wish to discuss some numerical checks we have made on these results. 

For our checks, we define an IR-regular differential cross section by sub- 
tracting the virtual infrared contribution as given by the YFS theory] 

IR " reg " 2a(ReE YFS ) J* , (4.18) 



dQkdkdQk dVLkdkdVL^ dVtkdkdVt 
with 5yfs defined as in Ref. [2^|. We have verified that our infrared limit 



A; — ^ is correct and in agreement with the YFS theory p0|| of that limit. 
This agreement with YFS theory is further illustrated in Fig. fj, where we 



show that ( [4.18|) is independent of our photon regulator mass 



m 



Further, looking into the results in Fig. 2, we see, in addition to their 
independence of lnm 7 , that the size of our 0(a 2 ) virtual correction is ~ 0.05 
of the respective real 0(a) correction itself, which is consistent with the 
naive LL power counting expectations, since ~ 0.05 here. Illustrations 
of the actual size of the subleading a 2 L part of the exact result ( p.3|) in 
the ALEPH SICAL-type acceptance in connection with BHLUMI4.xx will 
appear elsewhere 0. 

For completeness, we note that we have also evaluated the exact result 
for the s-channel exchange contribution to the 0(a 2 ) correction to the cross 
sections evaluated here. We find that it is below 0.04 of the 0(a 2 ) correction 
itself. This shows that such exchanges are indeed negligible at the level of 
precision of interest to us here in the evaluation of our respective 0(a 2 ) 
effects. 

Thus, the important exact 0{a) correction to e + e~ — > e + e~ +7 has now 
been calculated and it has been verified that it has the correct IR limit as 
well as a size consistent with the naive power counting expectations. The 
implementation of this exact result into the Monte Carlo event generator 
BHLUMI4.xx |j] for LABH is now in progress and will appear elsewhere 0. 
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Figure 2: lnm 7 dependence of the ratio R = ^Z^nJ dnkZl, where 
the numerator is defined in ( 4.18|) and the denominator is the 0{a) single 
bremsstrahlung cross- sect ion. We plot R as a function of k = 1 — s'/s and 
ln(m 7 /m e ) for positron polar angles (9, <p) = (2°, 0°) and photon polar angles 
(9,(f>) = (1.5°, 10°) with respect to the positron beam axis. This shows that 
lnm 7 has cancelled out when m 7 C m e , so that the result is independent of 
the infrared regulator. 
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5 Conclusion 



In this paper, we have presented the first complete exact 0(a) correction to 
e + e~ — > e + e~ + 7 at low angles in the SLC/LEP energy regime. We have 
checked the infrared limit of our work against known expectations. We have 
evaluated our results in the ALEPH SICAL type acceptance and found that 
the respective correction is consistent in size with the naive power counting 
expectations for its LL content. Implementation of our results into the Monte 
Carlo event generator BHLUMI M for LABH is in progress and will appear 
elsewhere 0. 

We should note that the authors in Ref. [ID have also given pioneering 
results for this 0(a) correction to e + e~ — > e + e~ + 7 at low angles. We have 
not compared our results with theirs however because they have not provided 
us with a final version of their fully differential results to date. We await their 
final published expressions with which we may compare. 

In summary, the dominant missing contribution to the 0(a 2 ) bremsstrah- 
lung effect in the theoretical prediction for the SLC/LEP luminosity process 
e + e~ — > e + e~, the 0(a 2 L) correction, has now been calculated exactly. We 
look forward with excitement to the application of this correction to reach 
a new level of precision on the respective theoretical prediction, via BH- 
LUMI4.xx H, of this luminosity. Such applications will appear elsewhereH. 
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Appendices 

A Scalar Integrals 

In this appendix, we define the scalar integrals used in the formulas in the 
text. We have, from Ref. [F7], with the kinematic conventions ( p. 14 ), 



B\2 


= B{m e 2 ;m 7 , m e ) 


(A.19) 


dM 
#13 


= B{m 2 — r;m 7 ,m e ) 


(A.20) 


-B24 


= B(t q ;m e ,m e ) 


(A.21) 


-B34 


= B(t p ;m e ,m e ) 


(A.22) 


r [r] 
^123 


= C(m 2 ,m 2 ,m 2 — r;m-y,m e ,m e ) 


(A.23) 


Cl24 


= Cim 2 ,t q ,m 2 \m^,m e ,m e ) 


(A.24) 


r [r] 


= C(m 2 - r,t p ,m 2 \m^,m e ,m e ) 


(A.25) 


C*234 


= C \m 2 ,tp,t q ] m e , m e , m e ) 


(A.26) 


D [r] 
-^1234 


= D(m 2 , m 2 , t p , m 2 , m 2 - r, t q ; m 7 , m e , m e , m e ) 


(A.27) 



where we have defined the basic scalar integrals as 

2 . (27T//) 4 -° f d D q 

B{p ,m l ,m 2 ) = / — , — — , 

7T z t J (q z — mf + te){{q + p) z — m 2 + te) 

(A.28) 

C(pl,pl, (pi +p 2 ) 2 ;m u m 2 ,m 3 ) = 

1 r d 4 q 
n 2 i J (q 2 — mi + ie)((q + p\) 2 — m 2 + ie)((q + p\ + p 2 ) 2 — m 2 + ie) 
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(A.29) 

d (Pi,P2,pI, (P1+P2+P3) 2 , (P1+P2) 2 , (P2 +P3) 2 ;m u m 2} m 3} m 4 ) = 



Tc 2 i J (q 2 — + ie) ((q + pi) 2 — m£ + ie) 



((<? + Pi +P2) 2 - rn£ + ze)((g + pi +p 2 + P3) 2 - W + «e) 

(A.30) 

The B integral is defined using dimensional regularization, and the C and 
D integrals are UV-finite. The final amplitude is independent of the mass 
scale jj, in the definition of B. 



B Mass Corrections 



In this appendix, we give our mass correction for the cross section in ( p. 17 ). 



Specifically, following Ref. JT5|, we find that the mass correction to (3.17) is 



dA(Tm ~ iP ° 2? '4n| 2 , (B.31) 



dVtkdkdVL k 2% 5 ss 
where we have defined 

& 2 Tfl 2 

\A m \ 2 = - j-^f (q-k, Pi ). (B.32) 

Here, the photon is radiated nearly parallel to q, and fo denotes the non- 
radiative cross section, summed over all polarizations, with the original q 
replaced by q — k. In the case of Bhabha scattering the Born cross section is 
proportional to the following invariant summed matrix element squared: 

f B e++e ~ = + (B.33) 

The complete non-radiative cross section for the 0(a 2 ) single bremsstrahlung 
mass corrections is then proportional to 

/o e++e " = (l + ^)f B e++e ~ (B.34) 
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with 



2 In 



- 1 



^B.35) 



From (p.32| ) it follows that, when summed over all fermion legs, the finite 
mass terms for the 0(a 2 ) single bremsstrahlung corrections are given by 



\A 



e++e~ |2 






2e e m 2 


1 + 


e 2 


{ qx ky 


4tt 2 


2e 6 m 2 


1 + 


e 2 


(q2k) 2 


4tt 2 


2e 6 m e 2 
(Pik) 2 


1 + 


e 2 
4vr 2 


2e 6 m e 2 


1 + 


e 2 


(p 2 k) 2 


4tt 2 



-T{-2 VlV2 + 2p 2 fc) 



(gigi ~ V\k f + (giPg ~ P2^) 2 

(<?1<?2 - <?2&) 2 

(P2?2 + P2k) 2 + {piq 2 + Pi A:) 2 

(<?i<?2 + gi/c) 2 
(pigi - qxk) 2 + (pigg - 92fc) 2 

(piP2 -p 2 k) 2 

(P2?2 + <?2^) 2 + (?lP2 + <?1^) 2 



(PlP2 2 



(B.36) 



This completes our appendices. 
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